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Abstract 

Nonlinear space-time dynamics, defined in terms of celebrated 'solitonic' equations, 
brings indispensable tools for understanding, prediction and control of complex behav- 
iors in both physical and life sciences. In this paper, we review sine-Gordon solitons, 
kinks and breathers as models of nonlinear excitations in complex systems in physics 
and in living cellular structures, both intra-cellular (DNA, protein folding and micro- 
tubules) and inter-cellular (neural impulses and muscular contractions). 
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1 Introduction 



In spatiotemporal dynamics of complex nonlinear systems (see [B El O S] ) , Sine-Gordon 
equation (SGE) is, together with Korteweg-deVries (KdV) and nonlinear Schrodinger 
(NLS) equations, one of the celebrated nonlinear-yet-integrable partial differential equations 
(PDEs)0 with a variety of traveling solitary waves as solution^ (see Figures [T] and [2l as 
well as the following basic references: O [H [71 [H [9l [TOl [11]). In complex physical systems, 
SGE solitons, kinks and breathers appear in various situations, including propagation of 
magnetic flux (fluxons) in long Josephson junctions [141115]. dislocations in crystals [16[ ll7j. 
nonlinear spin waves in superfluids [13] , and waves in ferromagnetic and anti- ferromagnetic 
materials [18\ [T9] - to mention just a few application areas. 

In this paper, we review physical theory of sine-Gordon solitons, kinks and breathers, 
as well as their essential dynamics of nonlinear excitations in living cellular structures, 
both intra-cellular (DNA, protein folding and microtubules) and inter-cellular (neural 
impulses and muscular contractions). 




-4-2 2 4 



Figure 1: Basic static examples of kinks: tanh(j;), arctan(3;) and bell-shaped soli- 
tons: sech(x), exp(— x^), together with their (absolute) differences; plotted in 
Mathematical . 

We show that sine-Gordon traveling waves can give us new insights even in such long- 
time established and Nobel-Prize winning living systems as the Watson-Crick double helix 
DNA model and the Hodgkin-Huxley neural conduction model. 

^ For a soft SGE-intro, see popular web-sites: [201 1211 [22] , Also, both KdV and NLS equations are 
mentioned in subsection 13.51 below as solitary models for muscular contractions on Poisson manifolds. 

■^A solitary wave is a traveling wave (with velocity v) of the form: 4>{x,t) — f(x — vt), for a smooth 
function / that decays rapidly at infinity; e.g., a nonhnear wave equation: (/>jj — 4>^^ = 4>i'^4>'^ ~ 1) ha-s a 
family of solitary-wave solutions: </>(2;,t) — sech(a;cosh/i + tsinh/^), parameterized by /i G R (see [H]). 
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2 Physical theory of sine— Gordon soUtons, kinks and breathers 



In this section, we give the basic theory of the sine-Gordon equation (and the variety of its 
travehng-wave solutions), as spatiotemporal models of nonlinear excitations in complex 
physical systems. 



2.1 Sine-Gordon equation (SGE) 

SGE is a real- valued, hyperbolic, nonlinear wave equation defined on R^'^, which appears 
in two equivalent forms (using standard indicial notation for partial derivatives: (f).^.^ = 
314) = d'^^jdx^): 

• In the (1+1) space-time (x, t)— coordinates, the SGE reads: 

= '/'xx - sin(/), or (/)tj(x,t) = (/>^^(x,t) - sin</>(x,t), (1) 

which shows that it is a nonlinear extension of the standard linear wave equation: 
4>tt — 4^ XX- The solutions 4){x,t) of ([T]) determine the internal Riemannian geometry 
of surfaces of constant negative scalar curvature R = —2, given by the line-element: 

ds2 = sin2(^) dt^ + cos^ dx^ 

where the angle (j) describes the embedding of the surface into Euclidean space 
(see [5]). A basic solution of the SGE ^ is: 



(j){x, t) = 4arctan 



, , X — vt 
exp ±- 



'VT 



(2) 



describing a soliton moving with velocity < v < 1 and changing the phase from 
to 27r {kink, the case of + sign) or from 27r to {anti-kink, the case of — sign)jl 

In the (1+1) light-cone (u, v)— coordinates, defined by: u = {x + t)/2, v = (x + t)/2, 
in which the line-element (depending on the angle between two asymptotic lines: 
u = const, V = const) is given by: 

ds'^ = dv? + 2 cos (jidudv + dv^ , 



^Each traveling soliton solution of the SGE has the corresponding surface in R'^ (see [12] 
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the SGE describes a family of pseudo-spherical surfaces with constant Gaussian 
curvature K = —1, and readsQ 

= sin (/>, or (j)^^{u,v) = sm(j){u,v). (3) 




Figure 2: Basic solitary SGE~solutions, simulated in Mathematica as systems of spring- 
coupled torsional pendula: (a) single soliton: 4){x,t) = 4 arctan ^exp ^Z^t ) ! (b) soliton- 

V sinh 



soliton collision: (j){x,t) = 4 arctan 



sinh —i^. 

4 arctan! ^^^^^^ 



; (c) soliton-antisoliton collision: (/>(x, t) = 

(modified 



cosh-ii^. 

1-v 

; and (d) single breather: — 4 arctan ( ^ qqsjJ 

and adapted from [IS])- 

A typical, spatially-symmetric, boundary-value problem for ([TJ is defined by: 

X £ [-L,L]CM, (teM+), 
,/.(x,0) = fix), Ux,0)=0, H-L,t) = ^iL,t), 

where f{x) G M is an axially-symmetric function (e.g., Gaussian or sech, see Figure [3]). 



''SGE ((21 is the single Codazzi-Mainardi compatibility equation between the first (Iq) and second (//c) 
fundamental forms of a surface, defined by the Gauss and Codazzi equations, respectively: 



Iq ~ du^ + 2 cos (f> du dv + dv'^ , 



He = 2 sin (j)dudv. 
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Figure 3: Numerical solution of the SGE ([T]) in Mathematica, using numerical ODE/PDE 
integrator NDSolve, with the following data (including the Gaussian initial state, zero 
initial velocity and symmetric boundary condition): x G [—20,20], t £ [0,30], (f){x,0) = 
exp(-a;2), ^i(x,0) = 0, (l){-20,t) = (/'(20,t). The waves oscillate around the zero plane 
and increase their width with time. Both near-periodicity and nonlinearity of the time 
evolution are apparent. 

Bdcklund transformations (BT) for the SGE ([T]) were devised in 1880s in Riemannian 
geometry of surfaces and are attributed to Bianchi and BackluncH (see, e.g. |25j). They 
have the form: 

-{4> + (p)^ = asm^—, -(0- v?)^ = -sm^— , (4) 

where both (f> and if are solutions of the SGE ([T|), and can be viewed as a transformation 
of the SGE into itself. BT (jH) allows one to find a 2-parameter family of solutions, given 

®In 1883, A. Backlund showed that if L : M — >■ M' is a pseudo-spherical line congruence between 
two surfaces M,M', then both M and M' are pseudo-spherical and L maps asymptotic lines on M to 
asymptotic lines on M' . Analytically, this is equivalent to the statement that if is a solution of the SGE 
m, then so are also the solutions of the ODE system (Q. 
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a particular solution 0g of ([T]). For example, consider the trivial solution cp = that, 
substituted into dH), gives: 



(f 2 . (f 



= — 2a sin — , f„ = sin 



a 



2' 



which, by integration, gives: 



2< = -21n(tan|) +p(r/), = -2 ln(tan ^) + 



from which the following new solution is generated: 

(/9 = 4 arctan 



exp(— rj + const) 

a 



The sine-forcing term in the SGE can be viewed as a nonlinear deformation: cj) — )• sin (/>, 
of the linear forcing term in the Klein-Gordon equation (KGE, a vacuum linearization of 
the SGE), which is commonly used for describing scalar fields (quantum) field theory: 



(5) 



This, in turn, implies that (as a field equation) SGE can be derived as an Euler-Lagrangian 
equation from the Lagrangian density: 

£sg(0) = ^(0?-</'')-1 + cos</). (6) 

It could be expected that £sg(</') is a 'deformation' of the KG Lagrangian: 

/:KG(0) = ^(0?-<^^)-f . (7) 

That can be demonstrated by the Taylor-series expansion of the cosine term: 



cos ( 



E 



(2n)! 

n=0 ^ ^ 



SO that we have the following relationship between the two Lagrangians: 
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Figure 4: Numerical solution of the SGE ()17p in Mathematical with the following data 
(including the Gaussian initial state, zero initial velocity and symmetric boundary condi- 
tion): x G [-20, 20], tE[0,30], (/)(x,0) =exp(-x2), 0) = 0, (/>(-20, t) = 0(20, t). 
Under the same boundary conditions, the SGE with the plus sine gives about 20 times 
higher amplitude waves, which are all above the zero plane and decrease their width with 
time. Again, both near-periodicity and nonlinearity of the time evolution are apparent. 

The corresponding Hamiltonian densities, of kinetic plus potential energy type, are 
given in terms of canonically-conjugated coordinate and momentum fields by: 



Both SGE and KGE are infinite-dimensional Hamiltonian systems [26|, with Poisson 
brackets given by: 




{F, G} = I 



oo 



6F 5G 5F 5G 



5(j){x) 5ir{x) 6'ir{x) 5(t){x) 



dx 



(8) 



— oo 
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so that both ([T]) and ([5]) follow from Hamilton's equations with Hamiltonian H and sym- 
plectic form a;|j 

= {H,(P}, 7rt = {H,7r}, 

n{(p,TT)dx, CO = dir Ad(j) dx. (9) 

-OO J ~oo 

The Hamiltonian Q is conserved by the flow of both SGE ([T]) and KGE (0), with an 
infinite number of commuting constants of motion (common level sets of these constants 
of motion are generically infinite-dimensional tori of maximal dimension). Both SGE and 
KGE admit their own infinite families of conserved functionals in involution with respect 
to their Poisson bracket ([8]). This fact allows them both to be solved with the inverse 
scattering transform (see 



2.2 Momentum and energy of SGE— solitons 

SGE is Lorentz-covariant (i.e., invariant with respect to special-relativistic Lorentz trans- 
formations; each SGE-soliton behaves as a relativistic object and contracts when v ^ c = 
the speed of light), and for this fact it has been used in (quantum) field theoryJll In 
Minkowski (1+1) space-time coordinates {x'^ G M^'^, x^ = t, x^ = x) with metric ten- 
sor r]^j^ (/^,z^ = 0,l;r]ii = — ?/22 = I'^ii = ^ii = 0)' the SG-Lagrangian density has 
the following 'massive form' of kinetic minus potential energy, with mass m and coupling 
constant A (see [^): 



1 




which reduces to the dimensionless form ([6]) by re-scaling the fields and coordinates: 



m 



(10) 



®For the Poisson-manifold generalization, see section |33] below. 
^The SGE, in both forms IT} and ([3]) has the following symmetries: 



t 


t + to, 


X — >• 


X, 


4>- 


(/!> (shift in t), 


t 




X X 


+ xo, 


<!>- 


(j> (shift in x), 


t 




X — ^ 


X, 




+ 2nn (discrete shifts in (p) 


t - 


-> -t, 


X — ^ 


X, 


<!>- 


-)■ (fi (reflection in t), 


t 




X — > 


-x, 


^- 


(j> (reflection in x), 


t 




X — > 


X, 




> —(j) (reflection in (f)), 


t 


t — vx 


X — > 


x — vx 




cj> (Lorentz transformations 









where e.g. reflection in <j) means: if <?!> is a solution then so is —0, etc. 
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The SG-Lagrangian density C^Q^^{(f>) = ^Csg{4>) obeys the conservation law and 
admits topological^ Noether current [with respect to ([TO]) ]: 

= — e^^du4> with zero-divergence: daj^ = 0, 
27r 

where e^'^ is the M^'^— Levi-Civita tensor. The corresponding topological Noether charge 
is given by: 

Q = j \dtf{x,t)\dx = ^\<t){+^,t)-<P{-^,t)\, 

with Qt = ^ |0i(-oo,t) - </>t(+oo,t)| = 0. 
zvr 

The most important physical quality of SGE is its energy-momentum (EM) tensor T^^,, 
which is the Noether current corresponding to spacetime-translation symmetry: — >■ 
x'^ + this conserved quantity is derived from the Lagrangian ([6]) as: 

T^,u = d/^cpducj) - ??^j,£sg(0)- 
T^^ has the following components [71 [3T]: 

Too = ^{(Pf + ^D + l- cos 0, Tio = = Toi. 
Til = ^(</.2 + 02)_i + cos0, 

EM's contravariant form T^*^ has the following components|£| 

TOO = Too, T^'=Tn, T^o = -Tqi. 

EM's conserved quantities are: momentum P = f T^^dx, which is the Noether charge with 
respect to space-translation symmetry, and energy E = J T^^dx, which is the Noether 
charge with respect to time-translation symmetry. Energy and momentum follow from 
EM's zero divergence: 

^rj.^.._^^j dtT^^ - 9,rio = ^{ dtE = dtjT^'dx = 
a^i - u ^ I ^^^01 _ Q^j^u ^0 ^\ dtP = dtjT^'^dx = ■ 



'The word topological means that it is not sensitive to local degrees of freedom, 
^The cont 
= 1/(^)m 



^The contravariant EM r'"^is obtained by raising the indices of T^^^ using the inverse metric tensor 
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2.3 SGE solutions and integrability 
2.3.1 SGE solitons, kinks and breathers 

The first 1-soliton solution of the SGE ([T]) was given by 

\/l — cos(a;t) 



in the form: 



{x,t) = 4arctan 



cosh(x\/l — w^) 



which, for w < 1, is periodic in time t and decays exponentially when moving away from 
X = 0. 

There is a weh- known travehng solitary wave solution with velocity v (see p^), given 
by the following generalization of ([2]): 



(x, t) = 4arctan 



exp ■ 



±2{z - zq) 



with {z = iJ,{x + vt)) , 



fill 



and the center at zq. In (jlip . the case +2 describes kink, while the case —2 corresponds 
to antikink. 

The stationary kink with the center at xq is defined by: 

4){x) = 2 arctan [exp(x — xq)] , 

(in which the position of the center xq can be varied continuously: — oo < xq < oo) and 
represents the solution of the first-order ODE: (p^ix) = sm(j){x). 

Regarding solutions of the slightly more general, three-parameter SGE: 

(f>tt = (^4>xx + bsm{X(t>), (12) 

the following cases were established in the literature (see [lOj and references therein): 

1. If a function w = (j){x, t) is a solution of ()12p . then so are also the following functions: 



27rn 



Wl 
W2 



±(f){Ci±x,C2±t) for (n = 0,±1,±2,...) 



±(/) ( X cosh C3 + t^/asmh. C3, x^^ — + t cosh C3 



where Ci, C2, and C3 are arbitrary constants. 
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2. Traveling- wave solutions: 



(x, t) = — arctan 
A 



exp lb 



if b\{Ci-aCl) > 0, 

(j){x,t) = - 



IT 4 

— + — arctan 
A A 



bXjCix + C2t + Cs) 
y/bXiCi - aCf) ^ 

bX{Cix + C2t + 



(13) 



exp lb 



y/bX{aC'l - C'i) 



if bX{Ci-aCf) < 0, 
where the first expression (for 6A(C| — aCf) > 0) represents another 1-soliton so- 
lution, which is kink in case of exp ( ^-^(pi^+^a^+Ca) j antikink in case of 

^ \^ ^bxicl-ac^) J 

_ b\{Cix+C2t+C3) \ ^^gg ^£ standard SGE (fTl), this kink-antikink ex- 

^bxicl-acf) J 

pression specializes to the Lorentz-invariant solution similar to pi|) : 



exp 



X,t) 



4 arctan 



exp 



±(X — Xr 



Vt 



(14) 



where the velocity f (0 < f < 1) and the soliton-center Xc are real-valued constants. 
The kink solution has the following physical (EM) characteristics: 

(i) Energy: 



(ii) Momentum: 



P[0K(x,t)] 

3. Functional separable solution: 

w{x, t) -- 



T'^dx 



T'^dx 



8v 



VT 



A 



arctan [f{x)g{t)] 



where the functions / = f{x) and g = g{t) are determined by the first-order au- 
tonomous separable ODEs: 



f^ = Af + Bf + C, 



9t 



-aCg^ + {aB + bX)g'^ - aA, 



where A, B, and C are arbitrary constants. In particular, for A = 0, B = > 0, 
and C > 0, we have the 2-soliton solution of 1281: 



w{x,t) = — arctan 
A 



7] sin(A;x -b Ai] 



ky/acosh{r]t + Bi 
where k, Ai, and Bi are arbitrary constants 



with (r/^ = ak"^ + bX > O) , 
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The only stable traveling wave SGE-solutions for a scalar field (p are 27r-kinks W\\ 
(localized solutions with identical boundary conditions = and (p = 27r). However, easier 
to follow experimentally are non-localized vr-kinks [43j (separating regions with different 
values of the field (/>), see also [18j and references therein. 

On the other hand, a breather is spatially localized, time periodic, oscillatory SGE- 
solution (see, e.g. [32]). It represents a field which is periodically oscillating in time and 
decays exponentially in space as the distance from the center x = is increased. This 
oscillatory solution of ([T]) is characterized by some phase that depends on the breather's 
evolution history. This could be, in particular, a bound state of vortex with an antivortex 
in a Josephson junction. In this case, breather may appear as a result of collision of a 
fluxon (a propagating magnetic flux-quantum) with an antifluxon, or even in the process 
of measurements of switching current characteristics, stationary breather solutions form 
one-parameter families of solutions. An example of a breather-solution of ([l]) is given by 

m- 

( sinT 

= 4 arctan 'n~~r^~^ 

\MCOsh [g[u)x) 

with parameters u = u{t) and T = T{t), such that 

rt 

g{u) = l/\/l + v? and T{t) = j g{u{t')) u{t') dt' . 

Jo 



2.3.2 Lax pair and general SGE integrability 

In both cases ([1]) and the SGE admits a Lax-pair formulation!^ 

L = [L,M], (15) 

where overdot means time derivative, L and M are linear differential operators and 
[L, M] = LM — ML is their commutator (or. Lie bracket). 



^'^The first Lax-pair for a nonlinear PDE was found by P. Lax in 1968 consisting of the following two 
operators [29) : 

L = -— —u, M = 4—7 — 6m- 3ux, 

ax-' ax-^ ax 



such that their Lax formulation (|15|) gives the KdV equation: 

ut — Quux + Uxxx — 0, by 

Lt = —Ut, LM — ML = Uxxx — Guux- 
The Lax-pair form of the KdV-PDE immediately shows that the eigenvalues of L are independent of t. 
The key importance of Lax's observation is that any PDE that can be cast into such a framework for other 
operators L and M, automatically obtains many of the features of the KdV-PDE, including an infinite 
number of local conservation laws. 
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Figure 5: Numerical solution of the damped, unforced SGE (|18p in Mathematica, with 
the following data (including the Gaussian initial state, zero initial velocity and symmet- 
ric boundary condition): x G [—20,20], t £ [0,30], (pix,0) = exp(— x^), (j)f(x,0) = 
0, (j){—20,t) = (/)(20,t), 7 = 0.2, F{x,t) = 0. Damping of the waves is apparent. 

For example, it was shown in [30J that the SGE ([T]) is integrable through the following 
Lax pair: 

0j = L(p, (j)^ = M(p, where (16) 



The Lax pair ()16p possesses the following complex-conjugate symmetry: if (j) = 
solves the Lax pair (|16p at (A,0), then ( ^ ) solves the Lax pair (|16p at {—X,(p). In 
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Figure 6: Numerical solution of the damped and spatially-forced SGE (fT8]l in 
Mathematica, with the following data (including the Gaussian initial state, zero ini- 
tial velocity and symmetric boundary condition): x G [—20,20], t G [0,30], (p{x,0) = 
exp(-x2), (/>t(x,0)=0, (/>(-20,t) = ^(20,t), 7 = 0.2, = 0.5 sech(x). We can see 

the central sech-forcing along all time axis. Damping of the SG-waves is also apparent. 

addition, there is a Darboux transformation for the Lax pair p6p as follows: let 



solves the Lax pair ()16p at (A,n). Also, from its spatial part: (p^ = Mcj), a complete 
Floquet theory can be developed. See [30] for the proofs and more technical details. 



ti = + 2i In 



(n G M) . 



If (j) = (j)\x=i, for some z/ G M, then 
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Figure 7: Numerical solution of the damped and temporally-forced SGE (jlSp in 
Mathematica, with the following data (including the Gaussian initial state, zero ini- 
tial velocity and symmetric boundary condition): x £ [—20,20], t G [0,30], </>(x,0) = 
exp(-x2), 0t(x,O)=O, (/>(-20,t) = 0(2O,t), 7 = 0.2, = 0.1 sin(t/2). We can see 

the sine-forcing along all time axis. Damping of the SG-waves is also apparent. 

2.4 SGE modifications 

2.4.1 SGE with the positive sine term 

The simplest SGE modification is to replace the minus sine term with the plus sine: 

4>tt = (t>xx + ^^T^'^(t>, or (|)^^{x,t) = (t)^^{x,t) +sm(t){x,t). (17) 

Again, a typical, spatially-symmetric, boundary-value problem for pT|) is defined by: 

X G [-L,L]CM, (tGM+), 
</.(x,0) = fix), Mx,0)=0, 4>{-L,t) = 4>{L,t), 

where f{x) G M is an axially-symmetric function (see Figure H]). 
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2.4.2 Perturbed SGE and vr kinks 



As we have seen above (and it was proved by [SIE]), the (1+1) SGE is integrable. In general 
though, the perturbations to this equation associated with the external forces and inhomo- 
geneities spoil its integrability and the equation can not be solved exactly. Nevertheless, 
if the influence of these perturbations is small, the solution can be found perturbatively 
[33j . The perturbation theory for solitons was developed by [57j and subsequently applied 
by [38] to dynamics of vortices in Josephson contacts. Perturbed SGEs come in a variety 
of forms (see, e.g. [371 [Ml EHl 112] ) . 




Figure 8: Numerical solution of the damped and both spatially and temporally forced SGE 
(I18p in Mathematica, with the following data (including the Gaussian initial state, zero 
initial velocity and symmetric boundary condition): x G [—20,20], t G [0,30], (j){x,0) = 
exp(-a;2), 4>t{x,0) = 0, (/'(-20,t) = (/)(20,t), 7 = 0.2, F{x,t) = 0.1sin(t/2) + 
0.5sech(j;). We can see both temporal sine-forcing and spatial sech-forcing along all time 
axis. Damping of the SG-waves is still visible. 

One common form is a damped and driven SGE: 

+ - (t^xx + sin = F, (18) 
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where ^(j)^ is the damping term and F{x, t) is the spatiotemporal forcing. Special cases of 
the forcing term F = F{x,t) in (fT8|) are: (i) purely temporal F = F{t) (e.g., periodic, see 
Figure [6]); (ii) purely spatial F = F{x) (e.g., central-symmetric, see Figure [7]); and (iii) 
spatiotemporal F = F{x,t) (e.g., temporally-periodic and spatially central-symmetric, see 
Figure ED. 

Considering (for simplicity) purely spatial forcing: F{x,t) = F{x), it has been shown 
in [m HH] that if F{xo) = for some point xq G M, this can be an equilibrium position 
for the soliton. If there is only one zero, in case of a soliton this is a stable equilibrium 

position if f^^j^) > 0; in case of an antisoliton, this is a stable equilibrium position if 



XO 



(^^aT^) ^ ^ ^^'^ references therein). 
In particular if 

F{x) = 2(/32 - l)sinh(^x)/cosh2(/3x), (;5 € M), 

the exact stationary kink-solution of (jlSp is: 

0;. = 4 arctan [exp {(3x)] . 

The stability analysis, which considers small amplitude oscillations around 
(pk [4>{k,x) = 0fc(x) + /(x)e'^*] , leads to the following eigenvalue problem: 

Lf = Tf, where L = -dl+[l-2 cosh'^ (/3x)] and F = - - 7A . 

The eigenvalues of the discrete spectrum are given by the formula 

F„ = /32(A + 2An-n2) - 1, 

where A(A + 1) = 2//3^. The integer part of A, yields the number of eigenvalues in the 
discrete spectrum, which correspond to the soliton modes (this includes the translational 
mode Fq, and the internal or shape modes F^ with n > (see [171 HH]). 

In case of a function F defined in such a way that it possesses many zeroes, maxima 
and minima, perturbed SGE psp describes an array of inhomogeneities. For example, 

9 „l3{x+x„) _ „3/3(x+x„) 

F(x)= V 4(1-/3^)" 



where Xn = (n + 2) log (\/2 + l) //3 (n = — —q + !••• ,(7 — l,g'), and g + 2 is the number 

lF(x 
dx 



of extrema points of F{x). When the soliton is moving over intervals where < 0, its 



internal mode can be excited. The points Xi where F{xi) = and '^^j^^^ < 0, are 'barriers' 
which the soliton can overcome due to its kinetic energy (for more details, see |47[ 
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Study of non-localized vr-kinks in parametrically forced SGE (PSGE): 



'/'it = (t>xx - <^{t/<i) sin</), (19) 

(over the fast time scale e, where a is a mean-zero periodic function with a unit amphtude), 
has been performed by [181 119j. via 27r-kinks as approximate solutions. In particular, a 
finite-dimensional counterpart of the phenomenon of vr-kinks in PSGE is the stabilization 
of the inverted Kapitza pendulum by periodic vibration of its suspension point. Geo- 
metrical averaging techniqu^ll] of [45j was applied as a series of canonical near-identical 
transformations via Arnold's normal form technique [51j, as follows. 
Starting with the Hamiltonian of PSGE (jl9p . given by: 

H{(f)) = J ( — + ~ a cos (f>j dx, where yp = (t>t = 

a series of canonical transformations was performed in [l8j with the aim to kill all rapidly- 
oscillating terms, the following slightly-perturbed Hamiltonian was obtained: 

which, after rescaling: X = ex, T = et, P = 2e~^p3, $ = 2^3, gave the following system 
of a slightly perturbed SGE with 27r-kinks as approximate solutions: 

<I>T = P + 0(e2), Pt = <l>xx - {all) sm + 0(e), 

where a_i is an anti-derivative with zero average. Finally, after rescaling back to variables 
(03, P3), approximate solutions <j)^ w Tp{x,t) in the form of vr-kinks were obtained, with 

'(/'(x, t) = 2 arctan ' 



exp eJ {ati 



where c is the wave-propagation velocity. For more technical details, see [18j . 

In addition, he following two versions of the perturbed SGE have been studied in |19j : 

1. Directly forced SGE: 

(t>tt - <t^xx + sin(/> = Mf{u)t). 
After shifting to the oscillating reference frame by the transformation: 

^ = e + Mio^^F{iot), (20) 



^^The averaged forces in a rapidly forced system (e.g. inverted Kapitza pendulum) are the constraint 
forces of an associated auxiliary non-holonomic system; the curvature of these constraints enters the 



expression for the averaged system [45] . 
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where F has zero mean and F"{t) = /{t), the parametrically forced ODE is ob- 
tained: 

e = - sin {9 + Muj-'^F{ujt)) , with (21) 

2 

H = ^ - A{cot)cos{9) + B{ujt)sm{e), 

where p is the momentum canonically conjugate to 9, and 

A{ujt) = cos{Muj''^F{ujt)), B{ujt) = sm( Muj~'^F{ujt)). 

Prom ()2ip . the corresponding evolution PDE (in canonical form) is obtained for a 
new phase 9 on top of a rapidly oscillating background field: 

9t =p, pt = 9^^ - sin {9 + Mu-'^F{ujt)). 

After retracing the identical transformation (j20p . the so-obtained (approximate) so- 
lutions become vr- kinks (see |19j for technical details). 

2. Damped and driven SGE 

0it - (l^xx + sin = Mfiujt) - acPt + r], (22) 



which is frequently used to describe long Josephson junctions [38] o Starting with 
a homogeneous transformation to the oscillating reference frame, analogous to ()20p 
and designed to remove the free oscillatory term: <j) = 9 + G{t), and substituting this 
transformation to (|22p . while choosing the function G so that it solves the ODE 



G + aG = Mf{ujt), 
the following evolution PDE is obtined (in canonical form) [19] : 

Ot = P, pt = 9^^-ap + 7]-sm{9 + G{ujt)). (23) 
For the particular case of /(t) = sinr, the the function G is found to be: 
^, , a M M . 

G[T) = K ^ cos T ^ Sm T. 

a; + + 

After a series of transformations (related to a directly- forced pendulum), in zeroth 
order in a, r], the evolution PDE (j23p reduces (after neglecting terms ~ w~^) to SGE, 
which has vr-kink solutions. Therefore, slightly perturbed vr-kinks are approximate 
solutions of the original equation ([22]) (on top of the rapidly oscillating background 
field; see [l9] for technical details). 



^^In (|22p. (f) represents the phase-difference between the quantum-mechanical wave functions of the two 
superconductors defining the Josephson junction, t is the normalized time measured relative to the inverse 
plasma frequency, x is space normalized to the Josephson penetration depth, while M f{u}t) represents 
tunneling of superconducting Cooper pairs (normalized to the critical current density). 
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2.4.3 SGE in (2+1) dimensions 

The (2+l)D SGE with additional spatial coordinate (y) is defined on M^'iasH 



Vtt = Av? - sin V9 = ip^^ + if - sin (p. 



(24) 



A special class of solutions of ()24p can be constructed by generalization of the solution of 
([T|) which does not depend on one of the coordinates, or, obtained by Lorentz transforming 
the solutions of a stationary 2D SGE. However, there are numerical solutions of (|24p which 
cannot be derived from the ([I]) or ([3]), e.g., radial breathers (or, pulsons) [34j- 

A more general class of solutions of the (2+l)D SG equation has the following formic 



which exactly satisfies (|24p with an arbitrary real- valued twice-differentiable function / = 
f{x lb t). The excitations, described by / are similar to elastic shear waves in solid 
mechanics |46| . 

Since the equation (j24p is Lorentz-covariant, we can obtain other solutions performing 
Lorentz transformations on ()25p , which leads to a class of solutions of the form |34] : 



""^^In the case of a long Josephson junction, the soliton solutions of (1241) describe Josephson vortices or 
fluxons. These excitations are associated with the distortion of a Josephson vortex line and their shapes 
can have an arbitrary profile, which is retained when propagating. In ((24]), (p denotes the superconducting 
phase difference across the Josephson junction; the coordinates x and y are normalized by the Josephson 
penetration length Aj, and the time t is normalized by the inverse Josephson plasma frequency (see 
|34) and references therein). 

^''Because of the arbitrariness of /, solution ()25p describes a variety of excitations of various shapes. 
Choosing / localized in a finite area, e.g., / = A/ cosh(a; — t), solution (|25|) describes an excitation, 
localized along x that keeps its shape when propagating, i.e., a solitary wave (in the sense of [?])• For 
each solitary wave of this type, there exists an anti-partner with an / of opposite sign in (|25p . For solitary 
waves to be solitons, there is an additional important criterion: restoring their shapes after they collide. 

Consider a trial function 



that, when t — oo, describes the propagation of two solitary shape waves toward each other (minus sign) 
or a solitary wave and its anti-partner (plus sign). One can see that (|25p can only approximately satisfy 
(|24|) when \ f'{x + t)f'{x — t)\ 1 for all values of x and t. This suggests that, in general, the condition 
for restoring the shapes may not be satisfied. In general case, (|24|) can not be satisfied, that prompts that 
the collision of two solitary waves leads to distortion of the original excitations |34j . 



^p{x, y,t) = A arctan exp [y — /(x ± t)] 



(25) 




'p{x, y,t) — 4: arctan exp [y — f{x + t) ± f{x — t)] 
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2.4.4 Two coupled SGEs 

The following two-parameter system of two coupled SGEs was introduced by 



-/32sin(</>-V), 



(26) 



ipii — a -0x2: = sm.{(f) — ■0), with constants (a,/3 > 0). 
For numerical solution, see Figure [9l 

10 



2x 10 




Figure 9: Numerical solution of the SGE-system ()26[) in Mathematical with the following 
data: xG [-10,10], te[0,5], a = 0.5, /3 = 0.3, (/)(x, 0) = 0.3exp(-x2), ijj{x,Q) = 
0.7exp(-x2), </>i(x,0)=0, (/>(-10,t) = 0(lO,t), Vt(2;,0)=0, ^(-10, t) = V'llO, t). 



The SGE-system (|26p has been exactly solved by j50J, where (using a series of substi- 
tutions) it was first reduced to the nonlinear second-order ODE: 



l + a2/32-c"(l + /3^ 



+ 



2lp 



12 



(c2 _ 1) (c2 _ a2) ^2 (^(e)2 + l 

equivalent to the following autonomous system in the {X, Y) —plane: 



(27) 



dX 



dY_ 



1 + a' 13' -c'{l + P'] 
(c2-l)(c2-a2)/i2' 



2y2 



X2 + 1 



X. 



(28) 
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System (|28p has an equihbrium point at the origin: (X, Y) = (0, 0), in which the Jacobian 
matrix is: 

. _/0 1\ , _ «^/3^-cMl + /3^) + l 

■^(O'") - U j ' ^ - (c2-l)(c2-a2)^2 • 

The phase portraits from this system show that there exist periodic solutions of the coupled 
SGEs (pUjl. 

Using the exponential ansatz: 

pexp(^) + gexp(-^) J. w ^m>\ 

<^{t,) = Tp—, r~?Y' ^i*'^ constants (p, q,r,s eR), 

rexp(^) + sexp(-4) 



four pairs of analytic solutions of the system (j27p . and therefore of the system (|26p . were 
found in [50j. We present here only the first two (simpler) solution pairs: 

Mx,t) = 77-^ sm(20 + ci^ + c2, 

= ——2 -T — 2 sin(2^) + ci^ + c2 - 2arctan(tan(^)), 

4 (c — i j /i 



^ = n{x-ct), c 



1-/32 ' 



/32 

02(x,t) = 771 — ^ sm(2^) + ci^ + c2, 

-02(2;, t) = — r^— — T-^sm(2^) + ci^ + c2-2arctan(cot(0), 

4 (_C ij /i 



^ = ^(x — ct), C 
For more technical details, see 1501. 



1-/32 • 



2.5 Sine— Gordon chain and discrete breathers 
2.5.1 Frenkel Kontorova model 

The original Prenkel-Kontorova model [W\ \W[ [T7j of stationary and moving crystal dis- 
locations, was formulated historically decades before the continuous SGE. It consists of a 
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chain of harmonically coupled atoms in a spatially periodic potential, governed by the set 
of differential-difference equations: 

'<Pn + [4>n+i - 2</'„ + 0„-i] + sin = 0, (29) 

where (p^ denotes the position of the nth atom in the chain. Alternatively, system (|29p 
represents a chain of torsionally-coupled pendula (see Figure [2]) , where is the angle 
which the nth pendulum makes with the vertical. 

2.5.2 Sine Gordon chain 

To derive dynamical equations of the sine-Gordon chain (SGC), consisting of anharmonic 
oscillators with the coupling constant fi, we start with the three-point, central, finite- 
difference approximation of the spatial derivative term cp^^ in the SGE: 

= --^[{(t>n-(l)n-l)-i<Pn+l-(l)n)]+0{x'^). 

Applying this finite-difference approximation to the SGE ([T]), and also performing the 
corresponding replacements: — )• i?!>„, (p^^ — )■ 0„ and ^ = we obtain the set of 

difference ODEs defining the SGC: 

(Pn + fJ' - <Pn~l) - (</'n+l " '^n)] + sin = . (30) 

The system (f30l) describes a chain of interacting particles subjected to a periodic on-site 
potential V{x) = sin(x). In the continuum limit, ([30]) becomes the standard SGE ([1]) and 
supports stable propagation of a kink-soliton of the form (jl4p . 




Figure 10: Simple sine-Gordon chain (SGC) with the coupling constant fi and the periodic 
on-site potential V{x) = sin(x). 

The linear-wave spectrum of (I30p around a kink has either one or two localized modes 
(which depends on the value of fi) [35J. The frequencies of these modes lie inside the 
spectrum gap. The linear spectrum, with the linear frequency lo and the wave number k, 
is given by: 

k 

uj^ = 1 + 4^810^ -, (31) 
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while the gap edge frequency is w = 1. 

The simplest example of (j30p . containing only two oscillators, is defined by |35j : 

''i>i,2 + fJ'{4'i,2 - '/'2,i) + sm</'i^2 = 0. (32) 

It was demonstrated in |44] . using the method of averaging in fast oscillations, that a 
perturbed SGC, damped and driven by a large-amplitude ac-force, might support localized 
kink solitons. Specifically, they considered the perturbed SGC: 

- [K+i - 20„ + K-i] + sin 0„ = X + a sin - 7</)„, (33) 

where x is a dc-force, a and u are the normalized (large) amplitude and frequency of a 
periodic force, respectively, while 7 is the normalized dissipative coefficient. Without the 
forcing on the right-hand side, (j33|) reduces to (pOj) . 



2.5.3 Continuum limits 



Perturbed SGEs have their corresponding perturbed SGCs. The following O-vr SGC was 
proposed in ^52] : 



+ 



sm 



(34) 



as an equation of a phase 0„-motion (of a O-vr array of Josephson junctions). Here, a is 
the lattice spacing parameter, 7 > is the applied bias current density, and 9n = (0 if 
n < and — vr if n > 0) is the phase jump of vr in The SGC equation ([3l|) is derived 
from the following discrete Lagrangian: 



E 



dt 



1 + COS (<?!>„ + 



dt. 



(35) 



In the continuum limit a <C 1 Lagrangian ()35p becomes 



1 + cos( 



dx dt . 



from which, the continuum limit of (j34p gives the following perturbed SGE: 



sin(0 + 61) + 7, 



where = (0 if x < and — vr if x > 0), while the differential operators Laipxx La(, 
are given by the following Taylor expansions: 



L. 



>n+l 



,2k 



k=0 



{2k + 2)1 



L. 



•-'n+l 



aH^x 



^ = E 



k=0 



{k+l) 



dzSina) 
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For more technical details, including several other continuum limits, see [52j . 
2.5.4 Discrete breathers 

Generally speaking, it is a well-known fact (see, e.g. [58] and references therein) that dif- 
ferent types of excitations, most notably phonons (propagating linear waves) and discrete 
breathers (DBs for short; they are time-periodic spatially localized excitations, also labeled 
intrinsic localized modes or discrete solitons) can occur as solutions of spatially-discrete 
nonlinear lattices. According to S. Flach et al. [Ml ES] DBs are caused by a specific 
interplay between the nonlinearity and discreteness of the lattice. The lattice nonlinearity 
provides with an amplitude-dependent tunability of oscillation or rotation frequencies of 
DBs, while its spatial discreteness leads to finite upper bounds of the frequency spectrum 
of small amplitude waves 1^ Although DBs present complex dynamical objects, experi- 
mental measurements can (in many cases) be well understood by using their time- averaged 
properties (see [56]). In addition, nonlinear discrete lattices admit different types of DBs 
depending on the spectrum of linear waves propagating in the lattice [531 ISH [551 157] , 
including: acoustic breathers, rotobreathers and optical breathers (see Figure [TT]). 




Figure 11: Different types of discrete breathers (DBs): acoustic breather (top), roto- 
breather (middle), and optical breather (bottom); modified and adapted from [531 [Ml [53 

m)- 

A particular system studied in |541 [59j has been characterized by the lattice Hamilto- 

^^DBs are not sensitive to specific types of nonfinearities in the iattice nor are tfiey confined to any 
lattice dimensions; tliey are (usually) dynamically and structurally stable and emerge in a variety of 
physical systems (ranging from lattice vibrations and magnetic excitations in crystals to light propagation 
in photonic structures and cold atom dynamics in periodic optical traps, see [59]). 
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nian: 

H = Y.^\xl + W{Xn-Xn-l) + V{Xn)) (36) 
n 

= Y.^\pl + W{Xn-Xn-l) + V{Xn)), 
n 

where x„ = are time-dependent coordinates with canonicahy-conjugate momenta 

Pn = Xn{t), W{xn) = W{x) is the nearest neighbor interaction, and V{xn) = V{x) is an 
optional on-site (substrate) potential. From (j36p the following equations of motion are 
derived: 

Xn = -W'{Xn- Xn-l) +W'{Xn+l- Xn) -V'{Xn), Or 

Xn = Xn, Pn = -W'{Xn - Xn-l) + W'{Xn+l - Xn) -V'{Xn), 

where (for simplicity) the following zero initial conditions are assumed: 

V{0) = W{0) = V'{0) = W'{0) = 0, V"{0) > 0, W"{0) > 0. 

Hamiltonian ()36p supports the excitation of small amplitude linear waves: 

Xn{t) ~ exp [i{ujgt - qn)] , 

with the wave number q and the corresponding frequency spectrum which, due to the 
underlying lattice, depends periodically on q: 

ul = V"{0)+4W"{0)sin'[l), 

and its absolute value has always a finite upper bound. The maximum (Debye) frequency 
of small amplitude waves is: 

cOg = ^/V"{0)+iW"{0). 

DBs exist for different types of potentials M^(a;) and V{x). DB solutions are [54] (i) 
time-periodic: x„(t + T^) = Xn{t), and (ii) spatially localized: x\n\^oo ~^ 0. In addition, 
if the Hamiltonian H is invariant under a finite translation/rotation of any — )• Xn, + A, 
then discrete rotobreathers may exist (see [60]), which are excitations characterized by one 
or several sites in the breather center evolving in a rotational state: XQ^t + Tf,) = xo(t) + A, 
while outside this center the lattice is governed again by time periodic spatially localized 
oscillations. 
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3 Sine— Gordon solitons, kinks and breathers in living cellular 
structures 



In this section, we give the apphcations of the sine-Gordon equation (and the variety of 
its travehng-wave solutions), as spatiotemporal models of nonlinear excitations in living 
cellular structures. 

3.1 SGE-solitons in DNA 

In this subsection, we review the first three papers describing SGE-solitary excitations in 
DNA0The first two papers in this domain were published by [62j and [64^1- incidentally, 
under the same title, in the same journal (PRA), using two slightly-different modifications 
of the same coupled SGE-system ([26]) . 

Firstly, Yomosa considered in [62] (see also [63]) the standard Watson-Crick double- 
helix i?— form DNA modelQ in which conformation and stability of DNA and the polynu- 
cleotide double helices are determined by0 

1. The energy Eb of the hydrogen H-bonds between inter-strand complementary base 
pairs, given by: 

EB = Y.B[l-COs{en-e'^-7T)], 

n 

where B is a parameter associated with the H-bond energy, while 9n = Z{Bn,Pn) 
and = ^{B'^, P^) denote angles between horizontal projections of the complemen- 
tary base pairs and their corresponding axes; and 

^''The idea that it is possible that soliton excitations may suggest a discovery of a new mechanism in the 
duphcation of DNA and the transcription of messenger ribonucleic acid (mRNA) goes back to [5T], who 
demonstrated the existence of transiently open states in DNA and synthetic polynucleotide double-helices, 
by hydrogen exchange measurements. 

^'^Note that, in the same period, Yakushevich et al. performed their SGE-solitary studies of DNA (see 
|66l 167) and references therein) , focusing on the effects of weak inhomogeneities in simple DNA fragments 
(consisting of uniform base sequences of a given type followed by uniform base sequence of the other type), 
which were described in terms of a parametrically-perturbed SGE. 

According to the Watson-Crick _B— form DNA model, the two polynucleotide strands forming a double 
helix are held together by hydrogen H-bonds. Yomosa was assuming that the H-bonding and the stacking 
energies (consisting of the electrostatic, the exchange, the charge-transfer, as well as the induction and 
dispersion interactions), were roughly proportional to the overlaps of molecular orbitals. 

^^Here, the zero-level of the energies Eb and Es are taken for the B— form of DNA and polynucleotide 
duplexes, while the mean energy of distorted double and triple H-bonds in ^ — T (adenine-thymine) and 
G — C (guanine-cytosine) base pairs is approximately represented (by the formula for molecular association 
in liquids by [65]) for the energy of a distorted single H-bond. 
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2. The stacking energy Es between intra-strand adjacent bases, given by: 

Es = S[l - cos{0n - On-i - ao)] + ^[l - cos(6'^ - 9'n_i - ao)], 



n 



where S" is a parameter associated with the stacking energy of DNA chains and 
ao = 36°. 

Next, by adding the rotational kinetic energy: 



n 



(where / is the mean value of the moments of inertia In of the bases for the rotations around 
the axes P) to the potential energies Eb and Es, the following SG-chain Hamiltonian for 
DNA and synthetic polynucleotide double-helices was formulated: 



n n 

+ J]{5[l-cos(0„-e„_i-Qo)] + 5[l-cos(0:,-Ci-ao)]}. (37) 



From the Hamiltonian (|37p . via canonical Hamiltonian formalism, the following two sets 
of coupled SGC-equations of motion were derived in [62j: 

I9n + B sin(6'„ - 6*^ - vr) + S'[sin(6'„ - 6n-i - ao) - sin(6'„+i - On - ao)] = 0, 
le'^ + B sin(0„ -e'^-Ti) + 5[sin(e; - - ao) - sH^'n+i - ^'n - "o)] = 0. 

Further, by linearizing this coupled ODE-system (assuming the smallness of the angles 
On — On-i — ao and 6'^ — 9'^_i — ao) and subsequently performing the continuum limit: 

Onit) ^ e{x,t), 9'^{t)^e'{x,t), (38) 

Yomosa derived the following system of two coupled SGEs, of the (slightly-modified) form 
of ([26]): 

I0tt-Se,, = -Bsm{e-e' -n), (39) 
le'u-Se'^^ = Bsm{9 -e' -tt). 

Unfortunately, in his time, Yomosa was not able solve the coupled system (I39p . so he 
took the difference of the two SGEs and obtained the following single SGE representing a 
dynamic (plane) base-rotator model: 

1 1 fs' I S 

■:^4>tt = 'Pxx-j2 " V 7' ^ " V 2B^' ^^^^ 
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Finally, by imposing the following boundary condition: 

cos</> = 1 for ((/) = 27rn, n = 0, ±1, ...) 
at ^ = iboo (at X = iboo for all t), 

the following traveling solitary-wave solutions of (|4U|) were obtained in the form of 
akink-antikink pair: 



4>{x,t) = 4 arctanexp 



Secondly, Zhang clarified in |64j the pioneering (and therefore somewhat-messy) ap- 
proach of Yomosa and proposed the following modified SGC-Hamiltonian for the form 
DNA double-helix (rewritten here in above Yomosa's notation for consistency): 

H = lY.^^(^n + On]+T.^{On,e'J (41) 

n n 

n 

where V{9n,0'n) is the inter-strand interaction energy in nth base pair, given by@ 

V{en,e'J = 5[l-cos(e„-e'„)] + A(l-cos0„) + A(l-cos0;) 
+ /3{3(l-cos0„cos0;)- [l-cos(0„-0;)]}. 

Prom the Hamiltonian (j4ip . the following SGC-equations of motion were derived in |64] : 

lOn + B sm{9n - 6l'„) + /3 [3 sin 6^ cos 9'^ - sm{9n - 9'J] + X sin 9n 

= S{9n+i-29n + 9n-i), (42) 

19'^ - Bsm{9n - 9'J + /3[3cos6l„ sin 9'^ + sin(6l„ - 9'J] + Asin6''„ 
= ^{9^+1 - 29'n + 9'n_i). 



By performing the approximation (|38p . Zhang introduced the continuum variables: 9 and 
9' . Subsequently, by introducing new variables: (p = 9 + 9' , tp = 9 — 9' , he obtained the 
following system of two perturbed and coupled SGEs: 

- i^/cmtt = (1/^') sin cP + {2/<f) sin(</./2) cos(V/2), 

- {l/cD^Pu = (1//2) sin V + { 2/<f ) sin (V>/2) cos((/)/2), 
where cq = ^/sjl, I = ^5/(3/3), d = y^S/X, 



■'Here the zero-level of the energy is taken for the B— form DNA, the same way as Yomosa did. 
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which cannot be solved analytically. So, by setting A = in (|12|) . Zhang arrived at the 
following system of two independent SGEs [(with Q = {2B + /3)/(3/3)]: 

ct>,,-{l/cl)4>,, = (l//2)sin</>, (43) 

V'xx - (l/co)V'tt = (Q//^)sinV, 

with the simple solution of a single soliton with velocity c: 

(/)^(a;, t) = 4 arctanexp (ibz) , ^/^q (x, t) = 4 arctan exp ^iL-^/Q , (44) 

where z = {-f - = x - ct, 7 = - c^/cq- 

Finally, by returning to original continuum variables 9 and 9' , from ()44p . Zhang obtained a 
set of vr-kink/antikink and 27r-kink/antikink solutions (see ^4] for more technical details). 

The third (and most-cited) paper in this domain (of SGE-solitary excitations in DNA) 
was |68) (see also |691I70)). who proposed a discrete SGC model for DNA-promoter dynam- 
ics. Salerno introduced the following SGC-Hamiltonian (slightly refined from Yomosa's 
and Zhang's Hamiltonians) : 

N N 

n=l n=l 
N 

ra=l 

where A'^ is the number of base-pairs in the SGC and K is the backbone spring constant 
along both DNA~helices. r]^ is a nonlinear parameter used to model the strength of H- 
bonds between complementary bases, chosen according to the following rule: 7/„ = An/3 
with A„ = 2 if it refers to ^ — T or T — ^ pairs, A^ = 3 otherwise, with (3 a free parameter. 
From the Hamiltonian (j45p . the following SGC-equations of motion were derived: 

^V'n = - 2'(/'„ + V-n-i) " sin(V'„ - 9n), (46) 

Wn = Ki9n+l-29n + 9n^l)- ^XnSm{9n-lpn)- 

Further, from ()46p . the following SGC-equation of motion is obtained for the angle differ- 
ence: (pj^ = i/j^ — 9n, between complementary bases: 

= '/'n+i - 20„ + 0„_i - — A„sin(/>„. (47) 

We remark that the ODE (j47p has the standard form of ()30p : from it, in the continuum 
limit, the standard SGE ([T|) is obtained, with exact soliton solutions (as described in 
the subsection 12.3.11 before). Salerno used the ODE-model ([17p to study nonlinear wave 
dynamics of the T7j41— DNA promoter (see |68j for further technical details). 
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3.2 SGE— solitons in protein folding 

For over a decade, it has been known that nonUnear excitations can influence conforma- 
tional dynamics of biopolymers; e.g., the effective bending rigidity of a biopolymer chain 
could lead to a buckling instability |115j . Subsequently, several models have been proposed 
to explain such protein transition (see, e.g. |116j and references therein). 

In this subsection, we review two protein-folding dynamics papers. Firstly, it was 
suggested in |114| that protein folding may be mediated via interaction of the protein 
(molecular) chain with SGE-solitons which propagate along the chain. Local potential 
energy of the chain is modeled by an asymmetric double-well potential: 

where the scalar variable (p represents local conformation of the protein, e is a small 
positive parameter, while 5 is the asymmetry parameter (ranging from —1 to 1). The two 
minima of the potential, corresponding to the local a- and /^-conformations of the chain, 
are positioned at tp = ±1 and the energy-difference between them is: /S.E = ^£(5. 

Solitonic excitations are realized in |114j by an additional, dissipative SGE-field ipix), 
where x is the position along the protein. The following interaction potential (with the 
positive parameter A) is used to mediate interaction between the two fields: 

The following dissipative equations of motion are derived [114j : 

(l>tt = (t>xx - sin(/. - 7^(/)t, (48) 

2A 

T^'Ptt = -4e(93 + '5)((/'^ - 1) - -cos(/)) -7<^V7i, 

where 700f and ^^p^Pt are dissipative terms. In the small interaction limit (ignoring 'y^(j)t), 
it is chosen: 

(f)(x, t) = f{x - vt) + Ae{x, t), ip{x, t) = 1 + A(p{x, t), 

where f{z) = 4arctan(e~^) is the usual SGE~kink, moving with velocity v. For ip{x,t), 
the following approximate sech-soliton solution is obtained: 

4A/m 1 



v^ + {uj/2f cosh^{x -vt)' 
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so, near the center of the (j) kink, ip is pushed away from its local minima 93 = 1 towards 
the other local minima. A localized static solution of (1481) is found to be: 



4 tan 



1 



^2 = 1 



A 



1 



e(l + A) 1 + ^72' 



where q = \/l — a sinh(2; — xq), a = h? /\^e{l + A)^]; for more technical details, see |114j . 

More recently, a Lagrangian field-theory based modeling approach to protein folding 
has been proposed in |117) . They proposed the protein Lagrangian including three terms: 
(i) nonlinear unfolding c/)'^— protein at the initial state: 



1 



I Aa, 



1 — COS 



(ii) nonlinear sources injected into the backbone, modeled by ip^ self-interaction: 

(iii) the interaction term (with the coupling constant A): 

Ant = -A(</'V)(V'^V')- 

The total potential (from all three terms) reads: 



Kot(</',V') 



A. 



1 — cos 



+ ^(V'W-A(<^V)(V'V). 



Assuming that X^pis small enough to be approximately at the same order with , the first 
term can be expanded in term of y'^A^j giving (up to the 2"*^ order accuracy): 

^tot(</', V') « ^ <AV - 4^ (<AV)' + ii^^i'f - A (</'V)(V'^V'), 

from which the total Lagrangian: £tot = ^uni + ^svc + ^int can be (up to the 2"'^-order 
accuracy) approximated by: 



£tot((A, ■0) 



2 

2 ^ ^ 4! ' ' 4! 



1 

2 



(49) 
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From the Euler-Lagrangian PDEs for the total Lagrangian ()49p : 

c^-^tot p. dCtot _ „ <9£tot o dCtot _ „ 

the following coupled and perturbed SGE and (nonlinear) KGE with cubic forcing are 
derived: 

sinf vi5|0|) +2A|,/,||V;|2, (50) 



^tt = V'..-^H' + 2A|V||</>P, (51) 

where A^— and A^— terms determine nonlinearities of backbone and source, respectively. 

Numerical solution of the two coupled (1+1) PDEs, ([50|) -([5T ]) . with the following 
boundary conditions: 

0(0, t) = (t){L, t) = and V(0, t) = '4){L, t) = for < t < 6, 

(j){x,^) = p{x) and ^{x,{)) = f {x) for < x < L, (52) 

4>^{x,Q) = q{x) and il)^{x,Q) = g{x) for < x < L, 

(where p(x), g(x), f{x) and g{x) are auxiliary functions), would describe the contour of 
conformational changes for protein folding. It was performed in [117] using the following 
forward finite differences: 



"^0 ■ I \/^4 
■ — — ^ sm 



where the values at the first time-step ti are given by the boundary conditions (I52p . while 
the values at the second time-step ^2 are determined by the 2"''^-order Taylor expansion: 



A% + -At ( + 2Af^pi 



sin ^^pi , (for i = 2,3,--- , - 1) 



V',^ = f^-Atg, + lAt'(l^±l^^±^ + 2ApU^-^f! 



For more technical details, see |117j . 
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3.3 SGE— solitons in microtubules 



In this subsection, we review two most-significant papers describing SGE-solitary exci- 
tations in microtubules (MTs)l^ and then point-out to some related quantum studies of 
neural MTs. 

3.3.1 Soliton dynamics in MTs 

To the best of our knowledge, the first paper describing soliton dynamics in MTs was 
[72j (see also [73l [TU [75] ; for a recent review, see ^j), in which Sataric et al. developed 
a ferroelectric model of neural MTs where the motion of MT subunits is reduced to a 
single longitudinal DOF per dimer at a position n, denoted by (p^. The overall effect of 
the surrounding dimers on a dipole at a chosen site n can be qualitatively described by 
the following double-well quartic potential: 

V{<pj = -^AcPl + \BcPt,, (53) 

where A and B are real parameters {A is a linear function of temperature) . As an electrical 
dipole, a dimer in the intrinsic electric field of the MT acquires the additional potential 
energy given by: 

Vei = -qE(f>n , 

where q = 18 x 2e {e = electron charge) denotes the effective mobile charge of a single 
dimer and E is the magnitude of the intrinsic electric field. In addition, two more MT- 
related energies can also be defined: 

(i) the potential energy of restoring strain-forces between adjacent dimers in the protofil- 
ament (with a unique spring/stiffness constant k): 

Vstr = ^H(t>n+i - <t>nf^ and 

^^MTs are major cytoskeletal proteins assembled from the tubulin protein that plays a crucial role in 
all eukaryotic cells. MT functions include cellular orientation, structure and guidance of membrane and 
cytoplasmic movement, which are crucial effects to cell division, morphogenesis, and embryo development. 
MT structure is a hollow cylinder that typically involves 13 protofilaments (of protein dimers called tubu- 
lins.). Each protofilament is formed from tubulin molecules arranged in a 'head-to-tail joint' fashion. The 
inner and the outer diameters of the cylinder are 15nm and 25nm, while its length may span dimensions 
from the order of micrometer to the order of millimeter. Each dimer is an electric dipole whose length 
and longitudinal component of the electric dipole moment are 8nm and 337Debye, respectively. The con- 
stituent parts of the dimers are a and /3 tubulins, corresponding to positively and negatively charged sides, 
respectively (see [7l]i as well as references in [72l[73l[74l[75] '). 
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(ii) the (velocity 0„— dependent) kinetic energy associated with the longitudinal displace- 
ments of constituent dimers with unique mass m: 

1 ^ 2 

n=l 

where is the total number of constituent dimers in the microtubular chain (MTC). The 
full MTC-Hamiltonian is now given by: 

H = T(0J + + ^'AJ + (54) 

Also, in order to derive a realistic equation of motion for the MTC, it is indispensable to 
include the viscous force: F^{(j)^) = — 7i?5>„, with the damping coefficient 7@ 

Using the Hamiltonian together with the damping force Fv ((/>„), and subsequently 
performing the continuum limit (with equilibrium spacing ro between adjacent dimers): 

4>n+i{t) (l){x,t) +ro(f>^{x,t) + -r'^(t)^^{x,t) + 

- Sataric et al. finally derived their nonlinear (forced and damped) wave equation: 

"i'/'ti + T-Z-t - krlcj)^^ = A(I)- Bcj}^ + qE. (55) 

The importance of the electric force term qE lies in the fact that the PDE (I55p admits 
soliton solutions with no energy loss, which acquires the form of a traveling wave, and can 
be expressed by defining a normalized displacement field [77| : 

where vq = y/k/mro is the sound velocity and v is the soliton-propagation velocity. In 
terms of the ■0(0 variable, the wave equation (j55p reduces to a damped anharmonic 
oscillator ODE: 

■0" + pV''^ + ■0 + o" = 0, where (56) 
p = -fv[m\A\{vl - v"^)]-^ , a = q\fB\A\-^/'^E, 

^■^We remark that the friction 7 should be viewed as an environmental effect, which does not lead to 
energy dissipation; this is a well-known result, relevant to energy transfer in biological systems [77] . 
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which has a unique bounded solution [72j: 



■0(0 = a H 7 \-) such that: (57) 

1 + exp (^e) 

_ a){i; -b){^-d)=^'-i;- (^f^^^ ■ (58) 

So the kink propagates along the protofilament axis with fixed velocity: 

V = vo/^l + 2-/^/{9d^m\A\), 



which depends on the strength of the electric field E via (j58p . The total conserved energy 
of the kink (I57p is given by 



|23 

2^/2^2 ^2 A 1 ^ 2 
E = — — — + —k— + -m V , where (59) 

6 ID 6 D 1 

4 niAa . , . , / a- ,. 
m = — ^ — IS kmk s etjective mass. 

For the further development of the theory, with kink-antikink waves traveling in opposite 
directions along the MTC, see [72l [78l [73l [H ES] . 

Now, from our general SGE perspective, Sataric model ([551) can be approximated by 
the perturbed SGE (|18p . rewritten here for our readers' convenience: 

(l^tt + l^l^t - (t^xx + sin = F, 

- if we apply the following assumptions: 

(i) normalized units: m = k = r^ = 1; 

(ii) electrical force: qE = F = F{x,t); and 

(iii) using the (first two terms of the) Taylor-series expansion of the sine term: 

sm4> = (I) --^^ + 0((f)^) ^ A4>- B4>^. (60) 
6 

Using assumptions (i)-(ii) and approximation ([60]) . all results of the section [2.4.21 (illus- 
trated with the simulation Figures HHS]) are ready to be employed for the further SGE 
analysis of solitary excitations in microtubules. 

The second paper describing soliton dynamics in MTs was ^79j . Under similar assump- 
tions as [72], they proposed the same expresions for the kinetic energy T{(l)^) and potential 



^■^The first term of (|59[1 expresses tlie binding energy of tlie kink and tlie second tlie resonant transfer 
energy (not tliat in realistic biological models, the sum of tfiese two terms dominate over tfie tliird term, 
being of order of leV). 
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energy V^tr of restoring strain-forces between adjacent dimers. However, in contrast to the 
quartic potential ([53]) . Chou et al. fohowed the recipes from soHd state physics [80J and 
expressed the interaction for the nth tubuUn molecule of a protofilament by the following 
periodic (effective) potential: 



V{cf>n) = Vo 



1 



COS 



ao 



where Vq is the half-height of the potential energy barrier, is the displacement of the nth 
tubulin molecule within a particular protofilament and ao = 8nm is the distance between 
the centers of two neighboring tubulin molecules along a protofilament. In this way, they 
defined the following MTC-Hamiltonian: 



H 



N 

E 

n=l 



(61) 



^2 + yofl-cos^ 
ao 



Using canonical Hamilton's equations, from (|6ip they derived the MTC~equations of mo- 
tion for the nth tubulin molecule within a protofilament (for n = 1,2, ...N): 



mcpn = k{4>n+l - ^(Pn + 4>n-l) Sm . 



ao 



ao 



(62) 



In the continuum limit: </>„(t) — s- </>(x,t), the MTC-equations of motion (j62|) reduce to the 
SGE: 



t«<44>^^ - ( Vosin^, or 



(63) 



1 



1 



r,2 p 



..,12 kal 1 47r2Vb 
sm0, with I c = — -r = , . , , 

ni kuq ' 



which is the same SGE as (I43|) that was used by [M] for DNA-solitons, with kink-antikink 
solutions (|44p . Using (j63p . Chou et al. showed that there was a very high and narrow peak 
at the center of the kink width, implying that a tubulin molecule would have its maximum 
momentum when it reaches the top of the periodic potential (for more technical details, 
see [791). 

In particular, in the case of neural MTs, possibility for sub-neuronal processing of 
information by cytoskeletal tubulin tails has been proposed by |81j . by showing that local 
electromagnetic field supports information that could be converted into specific protein 
tubulin-tail conformational states. Long-range collective coherent behavior of the tubulin 
tails could be modeled in the form of sine-Gordon kinks, antikinks or breathers that 
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propagate along the microtubule outer surface, and the tubulin-tail soliton collisions could 
serve as elementary computational gates that control cytoskeletal processes. The authors 
of [81] have used the results of [82], combined with the elastic ribbon SGE-model of 
|80j . Applying Backlund transformations (jH they found 2- and 3-soliton solutions, as 
well as their elastic collisions. They developed Map^e^^^-based animations of a whole 
'zoo' of colliding solitons, including kink/antikink pairs and three types of breathers: (i) 
a standing breather, (ii) a traveling large amplitude breather and (iii) a traveling small 
amplitude breather o 



3.3.2 Liouville's stringy time arrow in neural MTs 

Now, recall that the term cubic in %[) in the equation of motion (I56p was responsible for 
the appearance of a kink-like classical solution. Let us formulate a Liouville (1 + 1)— string 
theory of the neural MT-complex, following [78], and consider a general polynomial in T 
equation of motion for a static tachyon: 

T"{i)+pT'{i) = P{T), (64) 

(where ^ is some space- like coordinate and P{T) is a polynomial of degree n), in which the 
'friction' term T' expresses a Liouville derivative. In this interpretation of the Liouville 
field as a local scale on the world-sheet it is natural to assume that the single-derivative 
term expresses the non-critical string /3-function, and hence is itself a polynomial R of 
degree m: T'(^) = R{T). Such equations lead to kink-like solutions [83] : 

T(^) = -^{sgn(a2a4)a2tanh[^a2(2; - vt)] - 02}, (65) 

Z(24 Z 
A3 - 30204 

where v = is the velocity, 

04 

which is a universal behavior for biological systems [77], showing the existence of a scheme 
which admits kink-like solutions for energy transfer without dissipation in cells. The struc- 
ture of the equation (|64p . which leads to kink- like solutions (j65p . is generic for Liouville 
strings in non-trivial background space-times. 

According to the conventional Liouville theorem^^ written here in terms of Poisson 



^■^The standing breather soliton could be obtained in vivo in experiments in which the electric-field 
vector acts perpendicular to the microtubule z— axis. If a local vortex of the electromagnetic field is 
created somewhere in the neuron, then the exerted action of the electric field vector along the z-axis 
will be zero and no traveling soliton would be born. This standing breather, swinging at certain tubulin 
tail could catalyze attachment /detachment of microtubule associated proteins and promote or inhibit the 
kinesin walk. 

^^Liouville theorem (|66p is usually derived from the continuity equation 

Pt + div(px) = 0. 
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brackets {•,.}: 

p^ = -{p,H}, (66) 

the phase-space density of the field theory associated with the matter DOF of the MT- 
complex evolves with time as a consequence of phase-space volume-preserving symmetries. 
More generally, statistical description of the temporal evolution of the MT-complex using 
classical density matrices p{4>^,t) |84j: 

Pt = -{p,H} + P'Gijdp^p, (67) 

where pi are momenta canonically-conjugate to the fields (^*, and Gij is the metric in the 
space of fields {(/>*}. The non-Hamiltonian term in (j67|) leads to a violation of the Liouville 
theorem (j66p in the classical phase space {(p^ ,Pj^, and constitutes the basis for a dissipative 
quantum-mechanical description of the system [SI] , upon density-matrix quantization. In 
string theory, summation over world sheet surfaces will imply quantum fluctuations of the 
string target-space background fields 

Using Dirac's quantization rule: {., .} — > — i [., .], the quantum version of ()67p reads 
(in terms of the quantum commutator [.,.]) [84j : 

ft = i[p,^]+i/3*G,,[^\p], (68) 

where the hat denotes quantum operators, and appropriate quantum ordering (in the sense 
of [86]) is understood. In ()68p . H is the Hamiltonian evolution operator, while 

p = ^P(a)|M'a)(M/,|, with (Tr(^) = l), 

a 

is von Neumann's density matrix operator, in which each quantum state occurs with 
probability P{a); von Neumann's entropy is defined as: 

S = -Tr(p[lnp]). 

The very structure of the quantum Liouville equation (j68p implies the following prop- 
erties of the MT-complex: [HI [851 [87] : 

(i) conservation of probability P: 

Pt = Tvtip) = 0; 

(ii) conservation of average energy {E)f. = TitipE) = {piP^)t = 0; and 

(iii) monotonic increase in entropy S: 

St = -Txtiplnp) = S{(3'GijP^) > 0, 
- which naturally implies a microscopic arrow of time within the MT-complex [87| . 
■^^However, the quantum energy fluctuations: 

SE = [«H'^ » -{«H 
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3.4 SGE— solitons in neural impulse conduction 

Recently, two biophysicists from the Niels Bohr Institute in Copenhahen, T. Heimburg 
and A. Jackson (see [88\ I89j). challenged the half-a-century old electrical theory of neural 
impulse conduction, proposed by A.L. Hodgkin and A.F. Huxley in the form of their 
celebrated HH equations (1963 Nobel Prize in Physiology or Medicine). The HH model, 
which relies on ionic currents through ion channel proteins and the membrane capacitor, 
is the presently accepted textbook model for the nerve impulse conduction. 

For our readers' reference, here is a brief on the HH model, which (in its basic form) 
consists of four coupled nonlinear first-order ODEs, including the cable equation for the 
neural membrane potential V, together with m, h and n equations for the gating variables 
of Na and K channels and leakage (see \92\ I93j: for recent reviews, see |94^ I95j): 

CmV = -gN.m^HV -Vn^.)-gKn\V -VK)-gL{V -Vl)+I]''\ 

m = -{am + bm)m + am, h = -{uh + bh) h + ah, (69) 

re = — (o„ + 6„) re + a„, where 

a„ = 0.1(y + 40)/[l-e-(^+^o)/iO], 6^ = 4e-(^+65)/i8, 

a„ = 0.01(F + 55)/[l-e-(^+55)/iO], 6„ = 0.125 e-(^+65)/80^ 

a„ = 0.07e-(^+65)/20^ 6„ = 1/[1 + e-(^+35)/i0]^ 

Here the reversal potentials of Na, K channels and leakage are: T4^a = 50 mV, Vk = —77 
mV and Vi^ = —54.5 mV; the maximum values of corresponding conductivities are: (7Na = 
120 mS/cm^, = 36 mS/cm^ and 5l = 0.3 mS/cm^; the capacity of the membrane is: 
Cm = 1 /iF/cm^. The external, input current is given by: 

/f * = ffsyn(K - Vc) - t^n), (70) 

n 

which is induced by the pre-synaptic spike-train input applied to the neuron i, given by: 
Ui(t) = Va ^„ 6(t — tin)- In (|70l) . tin is the reth firing time of the spike-train inputs, gsyn and 
Vc denote the conductance and the reversal potential, respectively, of the synapse, Tg is 
the time constant relevant to the synapse conduction, and a{t) = {t/Tg) e~*/'^=0(t) (where 
Q{t) is the Heaviside function). 

In addition, Hodgkin and Huxley assumed that the total current is the sum of the trans- 
membrane current and the current along the axon and that a propagating solution exists 

are time-dependent and actually decrease with time |78j : 

dt{5Ef = -i « [l3\HWGj, »=« P'Gj.j^P' » 

= - « Q^P'GrjP' »= - « Q^dtC » < 0. 
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that fulfills a wave equation, so the simple cable ODE figuring in the basic HH model (j69p 
was further expanded into the following PDE for the propagating nerve impulse depending 
on the axon radius a: 

■^Vxx = CmVt + gxiV - Ek) + g^^{V - ENa), 

where Ri is the resistance of the cytosol within the nerve (see [89] for technical review). 

The HH model was originally proposed to account for the property of squid giant axons 
|921I93] and it has been generalized with modifications of ion conductances. More generally, 
the HH-type models have been widely adopted for a study on activities of transducer neu- 
rons such as motor and thalamus relay neurons, which transform the amplitude- modulated 
input to spike-train outputs. For attempts to relate the HH model (as well as its simplified 
form, Fitzhugh-Nagumo model (FHN]£0 [96ll97]) with propagation of solitons in neural 
cell membranes see [98] and references therein. 

However, the HH model fails to explain a number of features of the propagating nerve 
pulse, includmg the reversible release and reabsorption of hea@ and accompanying 
mechanical, fluorescence, and turbidity changes [88^ I89j: 

"The most striking feature of the isothermal and isentropic compression 
modulus is its significant undershoot and striking recovery. These features 
lead generically to the conclusions (i) that there is a minimum velocity of a 
soliton and (ii) that the soliton profiles are remarkably stable as a function of 
the soliton velocity. There is a maximum amplitude and a minimum velocity 
of the solitons that is close to the propagation velocity in myelinated nerves..." 

Earlier work of A.V. Hill [91] (another English Nobel laureate), on heat production in 
nerves (which was based on his previous work on heat production in contracted muscles 
|90j ) is actually reviewed in [93], where it is noted that the heat release and absorption 
response during the action potential is important 'but is not understood'. 

Based on thermodynamic relation between heat capacity and membrane area com- 
pressibility, Heimburg and Jackson considered in [88\ l89] a (1+1) hydrodynamic PDE for 
the dispersive sound propagation in a cylindrical membrane of a density-pulse, govern- 
ing the changes Ap^ (along the x-axis) of the lateral membrane density p^, defined by: 

^^Among several forms of the FHN-model, the simplest one (similar to the Van der Pol oscillator) is 
suggested by FitzHugh [96] : 

e|| = ex = X — — y, 
d7 = y = ix-y + b, 
where x is voltage (the fast variable), y is the slow recovery variable and 7, 6, e (0 < e 1) are parameters. 

^^Electrical currents through resistors generate heat, independent of the direction of the ion flux. The 
heat production in the HH-model should always be positive, while the heat dissipation should be related 
to the power of a circuit through the resistor, i.e. Q = P = V- I>Q (for each of the conducting objects 
in all phases of the action potential; see [881 189] 1. 
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Ap'^(x,t) = p^{x^t) — Pq , where = 4.035 • 10^^ is the equihbrium lateral area 

density in the fluid phase of the membrane slightly above the melting point. The related 
sound velocity c can be expanded into a power series (close to the lipid melting transition) 
as: 



c2 = c2+p(A/) + ^(A/) 



A\2 



+ ... 



(71) 



where cq = 176.6 m/s is the velocity of small amplitude sound, while p and q are parameters 
{p = -16.6 cg//,^ and q = 7<d.^cl/{p^f). 

In our standard notation, with cf){x,t) = Ap^(x,t), the dispersive wave equation of 
[551 [S9j can be written as: 

0it = c'^,, - ficP) . (72) 
Here, we need to make two remarks regarding the dispersive wave equation (|72p : 

1. If the compressibility is approximately constant and if Ap^ ^ p^, then the dispersive 
force /((/>) is zero and ([72]) reduces to the standard wave equation (depending only 
on the small amplitude sound Cq): 



^09 XX • 



2. If higher sound frequencies (resulting in higher propagation velocities as the isen- 
tropic compressibility is a decreasing function of frequency) become dominant, the 
dispersive forcing function /(</>) in ([72]) needs to be defined, or ad-hoc chosen [88] [89] 
to mimic the linear frequency-dependence of the sound velocity with a positive pa- 
rameter h as: /((^) = hcpr^rcxx ■ this case, the expansion ([TT]) needs to be explicitely 
included into PDE ([72]) . resulting in the equation governing dispersive sound prop- 
agation, which reads (in original notation of |88] 189]): 



A A 



d_ 

dx 



{cl+p{Ap' 



(73) 



Furthermore, by introducing the sound propagation velocity v, after the coordinate 
transformation: z = x — v ■ t, the dispersive PDE ([73p can be recast into a time- 
independent form, describing the shape of a propagating density excitation: 



Ap- 



d_ 

dz 



{ct+p{Ap^) + q{Ap 



A\2\ 



d_ 

dz 



A/ 



This 1-dimensional PDE has a localized (stationary) solution [99] : 



A/(z) 



p 



1+1 l + 2W^33^cosh(^^ji 



which is a sech-type soliton, a typical solution for KdV and NLS equations. 
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Now, without arguing either pro- or contra- Heimburg-Jackson theory of neural sound 
propagation, as an alternative to Hodgkin-Huxley electrical theory, we will simply accept 
the natural solitary explanation of the nerve impulse conduction, regardless of the physical 
medium that is carrying it (sound, or heat, or electrical, or smectic liquid crystal ^98], or 
possibly quantum- mechanical |78] ) . However, we are free to chose a different form for 
the dispersive force term f{<j)) in the perturbed wave equation (j72p . For example, instead 
of the Heimburg-Jackson's ad-hoc choice of the forth-derivative term, following [72\ and 
subsequent studies of neural MTs, we can choose a double-well quartic dispersive potential 
(1531) . which will, combined with the approximation (fGOl) . result in the standard SGE ([T]), 
generating analytical solutions of the traveling soliton, kink-antikink and breather form 
(as described in the subsection 12.3.11 before). 



3.5 Muscular— contraction solitons on Poisson manifolds 

For geometrical analysis of nonlinear PDEs, instead of using common symplectic struc- 
tures arising in ordinary Hamiltonian mechanics, the more appropriate approach is a Pois- 
son manifold (g*, {F,G}), in which g* is a chosen Lie algebra with a (it) Lie-Poisson 
bracket {F,G}±{fj,)) and carries an abstract Poisson evolution equation F = {F,H}. 
This approach is well-defined in both the finite- and the infinite-dimensional case (see 

[nniEiiniira). 

Recall that a Lie algebra consists of a vector (e.g., Banach) space g carrying a bilinear 
skew-symmetric operation [, ] : g x g — )• g, called the commutator or Lie bracket. This 
represents a pairingl^ rj] = ^rj — r/^ of elements ^ € g and satisfies Jacobi identity: 

[[C,r/],^] + [[77,H,e] + [[/U,e],??] =0. 

Let g be a (finite- or infinite-dimensional) Lie algebra and g* its dual Lie algebra, 

that is, the vector space paired with g via the inner product <, >: g* x g — )• R. If g is 

finite-dimensional, this pairing reduces to the usual action (interior product) of forms on 

vectors. The standard way of describing any finite-dimensional Lie algebra g is to provide 

its structural constants 7^^-, defined by = JijCk^ some basis ^j, (i = 1, . . . ,n) 

For any two smooth functions i^, G : g* — )• M, we define the (it) Lie-Poisson bracket 
isol 

{F,G}±(/.) = ±</., 



> ■ (75) 



^^Let El and E2 be Banach spaces. A continuous bilinear functional <, >: £1 x £'2 — > K is nondegenerate 
if < x,y > = implies a: = and y = for all x £ Ei and y £ E2. We say Ei and E2 are in duality if 
there is a nondegenerate bilinear functional <,>: Ei x £2 — >■ K- This functional is also referred to as an 
L'^ —pairing of Ei with £2- 

^°The (±) Lie-Poisson bracket (|75|) is a bilinear and skew-symmetric operation. It also satisfies the 
Jacobi identity: 

{{F, G}, ff}±(M) + {{G, H}, F}±{p) + {{H, F}, G}±(/i) = 
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Here n £ K*, [^■,^A is the Lie bracket in g and 6F/6fj,, 6G/6fi G g are the functional 
derivativeqf^ of F and G. 

Given a smooth Hamiltonian function : g* — )• M on the Poisson manifold (g*, {F, G}± (fJ.)), 
the time evolution of any smooth function F : g* — )• M is given by the abstract Poisson 
evolution equation: 

F = {F,H}. (76) 

Now, the basis of Davydov's molecular model of muscular contractior@ is oscillations of 
Amid I peptide groups with associated dipole electric momentum inside a spiral structure 
of myosin filament molecules (see |1UH I1U2^ I1U3| ). There is a simultaneous resonant 
interaction and strain interaction generating a collective interaction directed along the 
axis of the spiral. The resonance excitation jumping from one peptide group to another 
can be represented as an exciton, the local molecule strain caused by the static effect of 
excitation as a phonon and the resultant collective interaction as a soliton. 

Davydov's own model of muscular solitons was given by the nonlinear Schrddinger 
equation (NLS) [Ml M El 

iVi = -V.. + 2x|V'lV, (77) 

for — oo < X < +00. Here ^ = ip{x,t) is a smooth complex-valued wave function with 
initial condition 'ip{x,t)\t=o = ipi^) and x is a nonlinear parameter. In the linear limit 
ix = 0) (!77p becomes the ordinary Schrodinger equation for the wave -f/^-function of the 
free ID particle with mass m = 1/2 |102j . 

To put this muscular-contraction model into a rigorous geometrical settings, we can de- 
fine the infinite-dimensional phase-space manifold V = {(V', G S(R, C)}, where 5(M, C) 
is the Schwartz space of rapidly-decreasing complex- valued functions defined on R). We 
define also the algebra x('P) of observables on V consisting of real-analytic functional 



(thus confirming that g* is a Lie algebra), as well as the Leibniz rule: 

{FG,H}±(^i) = F{G,H}±{t^) + G{F,H}±{^). (74) 



''^For any two smooth functions F : g* — >■ R, we define the Frechet derivative D on the space L(g*,R) 
of all linear diffeomorphisms from g* to R as a map DF : g* — >■ L(g*,R); /i 1— >■ DF{^). Further, we define 
the functional derivative SF/Sfi £ g by 

SF 

DFU) ■ 5^1 ^< 5fi,— > 
on 

with arbitrary 'variations' S/j. £ g*. 

■^^For general overview of muscular contraction physiology and mechanics, see [1051 1100] . 

•^•^For a different (financial) application, with a variety of traveling- wave solutions of the NLS ((TTJ, 
including both sech- and tanh-solitons, solved in terms of Jacobi elliptic functions, see [106l I107j and 
references therein. 
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derivatives 5F /5%l), SF/dtp £ ^(M, C). The Hamiltonian function : P — t- M is given by: 



/ + 00 



and is equal to the total energy of the soliton, which is a conserved quantity for ([77|) (see, 

e.g. [iMiiim]). 

The Poisson bracket on xi'P) represents a direct generalization of the classical finite- 
dimensional Poisson bracket [104J: 

{F,G}+(^) = i/ I - ]dx. (78) 



— oo 



6ip 6ip dil) SiIj 



It manifestly exhibits skew-symmetry and satisfies Jacobi identity. The functional deriva- 
tives are given by: 

5F/6i; = -i{F,i;} and 6F/64) = i{F,i;}. 

Therefore the algebra of observables x('P) represents the Lie algebra and the Poisson 
bracket is the (+) Lie-Poisson bracket {F,G}^{ip). The nonlinear Schrodinger equation 
()77p for the solitary particle-wave is a Hamiltonian system on the Lie algebra x('P) relative 
to the (+) Lie-Poisson bracket {F, G}+(V') and Hamiltonian function H^tp). Therefore 
the Poisson manifold (x('P), {F,G}^{ip)) is defined and the abstract Poisson evolution 
equation (|76p . rewritten here as: ip = {ip,H}, which holds for any smooth function 
4> ■ XCP) is equivalent to (1771) . 

An alternative model of muscular soliton dynamics is provided by the Korteweg- 
deVries equation (KdV, see [JLQ4J)Ell 

ft - Qffx + /xxx = 0, (/, = d^f), (79) 

where x G M and / is a real-valued smooth function defined on R. This equation is related 
to the ordinary Schrodinger equation by the inverse scattering method (see [1081 [109]). 

Again, we may define the infinite-dimensional phase-space manifold V = {/ S S'(M)}, 
where S'(M) is the Schwartz space of rapidly-decreasing real-valued functions M). Further, 
we define x(V) to be the algebra of observables consisting of functional derivatives 6F/6f G 
The Hamiltonian : V — t- M is given by: 

f+OO 



/+00 1 
^ if' + ^f!)dx, 



■^■^The most common traveling-wave solutions of the KdV (|79p are sech-solitons with the velocity c, of 
the form [12]: 



f{x, t) = -sech 
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and provides the total energy of the soHton, which is a conserved quantity for ()79p (see 



As a real-valued analogue to ([78j) . the (+) Lie-Poisson bracket on x(V) is given via 



Again it possesses skew-symmetry and satisfies Jacobi identity. The functional derivatives 
are given by 5F/5f = {F,f}. The KdV equation ([79]) . describing the behavior of the 
muscular molecular soliton, is a Hamiltonian system on the Lie algebra x(V) relative to 
the (+) Lie-Poisson bracket {F,G}+{f) and the Hamiltonian function H{f). Therefore, 
the Poisson manifold (x(V), {F, is defined and the abstract Poisson evolution 

equation (|76p . rewritten here as: / = {/,-?/}, which holds for any smooth function 
/ : X(V) -^M, is equivalent to (|79]l . 

Another alternative model of muscular soliton dynamics is provided by our SGE ([T|): 
(pn = (j)^^ — sin(/). Again, we may define the infinite-dimensional phase-space manifold 
V = {(/> G S'(R)}, where 5(R) is the Schwartz space of rapidly-decreasing real-valued 
functions M). Further, we define x(V) to be the algebra of observables consisting of 
functional derivatives 5F/d(f) G S'(M). The Hamiltonian : V — )• M is given by: 



and provides the total energy of the soliton, which is a conserved quantity for the SGE 
([1]). The (+) Lie-Poisson bracket on x(V) is given via ([71]) by: 



Again it possesses skew-symmetry and satisfies Jacobi identity. The functional derivatives 
are given by: 6F/6(f) = {F,(j)} G S'(M). The SGE ([1]), describing the behavior of the 
molecular muscular soliton, is a Hamiltonian system on the Lie algebra xO^) relative to 
the (+) Lie-Poisson bracket {F, G}+((/>) and the Hamiltonian function H((p). Therefore, 
the Poisson manifold {xiV),{F,G}+{4>)) is defined and the abstract Poisson evolution 
equation (j76p . rewritten here as: = {(j),H}, which holds for any smooth function 
</• '■ — ^K, is equivalent to ([I]). 

4 Conclusion 

In this paper, we have reviewed sine-Gordon equation and its traveling wave solutions. 
These solitary spatiotemporal processes can serve as realistic models of nonlinear excita- 
tions in complex systems in physical sciences as well as in various living cellular structures, 



[T08l[T09] l. 



m by: 





- (vr^ + (pi) + 1 — cos (f) dx, 
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including intra-cellular ones (DNA, protein folding and microtubules) and inter-cellular 
ones (neural impulses and muscular contractions). We have showed that sine-Gordon 
solitons, kinks and breathers can give us new insights even in such long-time established 
and Nobel-Prize winning living systems as the Watson-Crick double helix DNA model 
and the Hodgkin-Huxley neural conduction model. 
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